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Abstract 

Based on first principle calculations, we show that a family of nonmagnetic materials 
including TaAs, TaP, NbAs and NbP are Weyl semimetal (WSM) without inversion 
center. We find twelve pairs of Weyl points in the whole Brillouin zone (BZ) for each of 
them. In the absence of spin-orbit coupling (SOC), band inversions in mirror invariant 
planes lead to gapless nodal rings in the energy-momentum dispersion. The strong 
SOC in these materials then opens full gaps in the mirror planes, generating nonzero 
mirror Chern numbers and Weyl points off the mirror planes. The resulting surface 
state Fermi arc structures on both (001) and (100) surfaces are also obtained and 
show interesting shapes, pointing to fascinating playgrounds for future experimental 
studies. 
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Most topological invariants in condensed matter non-interacting phases are defined on 
closed manifolds in momentum space. For gapped systems, both the Chern insulator and 
Z 2 topological insulator phases can be defined using the Berry phase and curvature in either 
the entire or half of the two-dimensional (2D) Brillouin zone (BZ) respectively. [1, 2] A 
similar idea can be generalized to gapless metallic systems. In three dimensional (3D) 
systems, besides the BZ, an important closed manifold in momentum space is a 2D Fermi 
surface (FS). Topological metals can be defined by Chern numbers of the single-particle wave 
functions at the Fermi surface energies. [3-5] Such nonzero FS Chern number appears when 
the FS encloses a band crossing point — Weyl point — which can be viewed as a singular 
point of Berry curvature or “magnetic monopole” in momentum space. [6] Material with 
such Weyl points near the Fermi level are called Weyl semimetals(WSM). [7-9] 

Weyl points can only appear when the spin doublet degeneracy of the bands is removed 
by breaking either time reversal T or spacial inversion symmetry P (in fact, Weyl points 
exist if the system does not respect T ■ P). In these cases, the low energy single particle 
Hamiltonian around a Weyl point can be written as a 2 x 2 “Weyl equation”, which is half 
of the Dirac equation in three dimension. According to the “no-go Theorem” ,[10, ll]for 
any lattice model the Weyl points always appear in pairs of opposite chirality or monopole 
charge. The conservation of chirality is one of the many ways to understand the topological 
stability of the WSM against any perturbation that preserve translational symmetry: the 
only way to annihilate a pair of Weyl points with opposite chirality is to move them to 
the same point in BZ. Since generically the Weyls can sit far away from each other in the 
BZ, this requires large changes of Hamiltonian parameters, and the WSM is stable. The 
existence of Weyl points near the Fermi level will lead to several unique physical properties, 
including the appearance of discontinuous Fermi surfaces (Fermi arcs) on surface, [6, 8, 9] 
the Adler-Bell-Jackiw anomaly [7, 12-14] and others. [15, 16] 

The first proposal to realize WSM in condensed matter materials was suggested in Ref. 
8 for Rn 2 lr 207 Pyrochlore with all-in/all-out magnetic structure, where twenty four pairs of 
Weyl points emerge as the system undergoes the magnetic ordering transition. A relatively 
simpler system HgCr 2 Se 4 [6] has then been proposed by some of the present authors, where a 
pair of double-Weyl points with quadratic touching appear when the system is in a ferromag¬ 
netic phase. Another proposal involves a fine-tuned multilayer structure of normal insulator 
and magnetically doped topological insulators. [17] These proposed WSM systems involve 
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magnetic materials, where the spin degeneracy of the bands is removed by breaking time 
reversal symmetry. As mentioned, the WSM can be also generated by breaking the spatial 
inversion symmetry only, a method which has the following advantages. First, compared 
with magnetic materials, nonmagnetic WSM are much more easily studied experimentally 
using angle resolved photo emission spectroscopy (ARPES) as alignment of magnetic do¬ 
mains is no longer required. Second, without the spin exchange field, the unique structure 
of Berry curvature leads to very unusual transport properties under strong magnetic field, 
unspoiled by the magnetism of the sample. 

Currently, there are three proposals for WSM generated by inversion symmetry break¬ 
ing. One is a super-lattice system formed by alternatively stacking normal and topological 
insulators. [18] The second involves Tellurium or Selenium crystals under pressure [19] and 
the third one is the solid solutions LaBii_a;Sba;Te 3 and LuBii_a;Sba;Te 3 [20] tuned around the 
topological transition points. [21] In the present study, we predict that TaAs,TaP, NbAs and 
NbP single crystals are natural WSM and each of them possesses a total of 12 pairs of Weyl 
points. Compared with the existing proposals, this family of materials are completely stoi¬ 
chiometric, and therefore, are easier to grow and measure. Unlike in the case of pyrochlore 
iridates and HgCr 2 Se 4 , where inversion is still a good symmetry and the appearance of Weyl 
points can be immediately inferred from the product of the parities at all the time reversal 
invariant momenta (TRIM), [22-24] in the TaAs family parity is no longer a good quan¬ 
tum number. However, the appearance of Weyl points can still be inferred by analyzing 
the mirror Chern numbers (MCN) [25, 26] and indices [22, 27] for the four mirror and 
time reversal invariant planes in the BZ. Similar with many other topological materials, the 
WSM phase in this family is also induced by a type of band inversion phenomena, which, 
in the absence of spin-orbit coupling (SOC), leads to nodal rings in the mirror plane. Once 
the SOC is turned on, each nodal ring will be gapped with the exception of three pairs of 
Weyl points leading to fascinating physical properties which include complicated Fermi arc 
structures on the surfaces. 


CRYSTAL STRUCTURE AND BAND STRUCTURE 

As all four mentioned materials share very similar band structures, in the rest of the 
paper, we will choose TaAs as the representative material to introduce the electronic struc- 
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tures of the whole family. The experimental crystal structure of TaAs [28] is shown in 
Fig. 1(a). It crystalizes in body-centered-tetragonal structure with nonsymmorphic space 
group lAimd (No. 109), which lacks inversion symmetry. The measured lattice constants 
are 0=6=3.4348 A and c=11.641 A. Both Ta and As are at 4o Wyckoff position (0, 0, u) 
with 0=0 and 0.417 for Ta and As, respectively. We have employed the software package 
OpenMX [29] for the hrst-principles calculation. It is based on norm-conserving pseudo¬ 
potential and pseudo-atomic localized basis functions. The choice of pseudopotentials, 
pseudo atomic orbital basis sets (Ta9.0-s2p2d2fl and As9.0-s2p2dl) and the sampling of 
BZ with 10 X 10 X 10-grid have been carefully checked. The exchange-correlation func¬ 
tional within generalized gradient approximation (GGA) parameterized by Perdew, Burke, 
and Ernzerhof has been used. [30] After full structural relaxation, we obtain the lattice 
constants a=6=3.4824 A, c=11.8038 A and optimized u=0.4176 for As site, in very good 
agreement with the experimental values. To calculate the topological invariant such as MCN 
and surface states of TaAs, we have generated atomic-like Wannier functions for Ta 5d and 
As 4p orbitals using the scheme described in Ref. 31 

We hrst obtain the band structure of TaAs without SOG by GGA and and plot it along 
the high symmetry directions in Fig. 1(c). We find clear band inversion and multiple band 
crossing features near the Fermi level along the ZN, ZS and ES lines. The space group of 
the TaAs family contains two mirror planes, namely M^, My (shaded planes in Fig. 1(b) ) 
and two glide mirror planes, namely M^y, (illustrated by the dashed lines in Fig. 1(b)). 

The plane spanned by Z, N and F points is invariant under mirror My and the energy bands 
within the plane can be labeled by mirror eigenvalues ±1. Further symmetry analysis shows 
that the two bands that cross along the Z to N line belong to opposite mirror eigenvalues and 
hence the crossing between them is protected by mirror symmetry. Similar band crossings 
can also be found along other high symmetry lines in the ZNF plane, i.e. the ZS and NS 
lines. Altogether these band crossing points form a “nodal ring” in the ZNF plane as shown 
in Fig.2(b). Unlike for the situation in ZNF plane, in the two glide mirror planes {M^y and 
M-xy), the band structure is fully gaped with the minimum gap of roughly 0.5 eV. 

The analysis of orbital character shows that the bands near the Fermi energy are mainly 
formed by Ta 5d orbitals, which have large SOG. Including SOG in the first principle cal¬ 
culation leads to a dramatic change of the band structure near Fermi level, as plotted in 
Fig. 1(d). At first glance, it seems that the previous band crossings in the ZNF plane are 
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FIG. 1. Crystal structure and Brillouin zone (BZ). (a) The crystal symmetry of TaAs; (b) 
The bulk BZ and the projected surface BZ for both (001) and (100) snrfaces; (c) The band strnctnre 
of TaAs calculated by GGA without including the spin-orbit coupling; (d)The band structure of 
TaAs calculated by GGA with the spin-orbit coupling. 

all gaped with the exception of one point along ZN line. Detailed symmetry analysis reveals 
that the bands “2” and “3” in Fig. 1(d) belong to opposite mirror eigenvalues, indicating the 
almost touching point along the ZN line is completely accidental. In fact there is a small 
gap of roughly 3 meV between bands “2” and “3” as illustrated by the inset of Fig. 1(d). 
The ZNF plane then becomes fully gapped once SOC is turned on. 


TOPOLOGICAL INVARIANTS FOR MIRROR PLANE AND WEYL POINTS 

Since the material has no inversion center, the usual parity condition [22-24] can not be 
applied to predict the existence of WSM. We then resort to another strategy. As previously 
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mentioned, the space group of the material provides two mirror planes and My)^ where 
the MCN can be defined. If a full gap exists for the entire BZ, the MCN would directly 
reveal whether this system is a topological crystalline insulator or not. Interestingly, as 
shown below, if the system is not fully gaped we can still use the MCN to find out whether 
the material hosts Weyl points in the BZ or not. Besides the two mirror planes, we have two 
additional glide mirror plane {M^y and M^^y). Although the MCN is not well defined for the 
glide mirror planes, the Z 2 index is still well defined here as these planes are time reversal 
invariant. We then apply the Wilson loop method to calculate the MCNs for the two mirror 
planes and Z 2 indices for the two glide mirror planes. For a detailed explanation of the 
method, please refer to Ref. 5 and 32. The results are plotted in Fig.2(d), which shows that 
MCN is one for ZNF plane (My) and Z 2 index is even or trivial for ZXF plane(Mj;j^). Then, 
if we consider the (001) surface, which is invariant under the My mirror. The non-trivial 
helical surface modes will appear due to the non-zero MCN in ZNF plane, which generates 
a single pair of FS cuts along the projective fine of the ZNF plane (the a:-axis in Fig. 2(c)). 
Whether these Fermi cuts will eventually form a single closed Fermi circle or not depends 
on the Z 2 index for the two glide mirror plane, which are projected to the dashed blue lines 
in Fig. 2(c). Since the Z 2 indices for the glide mirror planes are trivial, as confirmed by our 
Wilson loop calculation plotted in Fig. 2(d), there are no protected helical edge modes along 
the projective lines of the glide mirror planes (dashed blue fines in Fig. 2(c)) and the Fermi 
cuts along the a:-axis in Fig. 2(c) must end somewhere between the a:-axis and the diagonal 
fines (dashed blue lines in Fig. 2(c)). In other words they must be Fermi arcs, indicating the 
existence of Weyl points in the bulk band structure of TaAs. 

From the above analysis of the MCN and Z 2 index of several high-symmetry planes, 
we can conclude that Weyl points exist in the TaAs band structure. We now determine 
the total number of the Weyls and their exact positions. This is a hard task, as the Weyl 
points are located at generic k points without any little-group symmetry. For this purpose, 
we calculate the integral of the Berry curvature on closed surface in k-space, which equals 
the total chirality of the Weyl points enclosed by the given surface. Due to the four fold 
rotational symmetry and mirror planes that characterize TaAs, we only need to search for 
the Weyl points within the reduced BZ - one eighth of the whole BZ. We first calculate 
the total chirality or monopole charge enclosed in the reduced BZ. The result is one, which 
guarantees the existence of, and odd number of Weyl points. To determine precisely the 
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FIG. 2. Nodal rings and Weyl points distribution, as well as Z 2 and MCN for mirror 
planes, (a) 3D view of the nodal rings (in the absence of SOC) and Weyl points (with SOC) in 
the BZ; (b) Side view from [100] and (c) top view from [001] directions for the nodal rings and 
Weyl points. Once the SOC is turned on, the nodal rings are gapped and give rise to Weyl points 
off the mirror planes (see movie in Appendix); (d) top panel: Flow chart of the average position 
of the Wannier centers obtained by Wilson loop calculation for bands with mirror eigenvalue i in 
the mirror plane ZNF; bottom panel: The flow chart of the Wannier centers obtained by Wilson 
loop calculation for bands in the glide mirror plane ZXF. There is no crossing along the reference 
line (the dashed line) indicating the Z 2 index is even. 

location of each Weyl point, we divide the reduced BZ into a very dense fc-point mesh and 
compute the Berry curvature or the ^^magnetic field in momentum space” [31, 33] on that 
mesh as shown in Fig. 3. From this, we can easily identify the precise position of the Weyl 
points by searching for the “source” and “drain” points of the “magnetic field”. The Weyl 
points in TaAs are illustrated in Fig. 2(a), where we find 12 pairs of Weyl points in the 
vicinity of what used to be, in the SOC-free case, the nodal rings on two of the mirror 
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TABLE I. The two nonequivalent Weyl points in the xyz coordinates shown in Fig. 1(b). The 
position is given in unit of the length of T-E for x and y and of the length of T-Z for z. 



Weyl Node 1 

Weyl Node 2 

TaAs 

(0.949, 0.014, 0.0) 

(0.520, 0.037, 0.592) 

TaP 

(0.955, 0.025, 0.0) 

(0.499, 0.045, 0.578) 

NbAs 

(0.894, 0.007, 0.0) 

(0.510, 0.011, 0.593) 

NbP 

(0.914, 0.006, 0.0) 

(0.494, 0.010, 0.579) 


invariant planes. For each of the mirror invariant planes, after turning on SOC, the nodal 
rings will be fully gaped within the plane but isolated gapless nodes slightly off plane appear 
as illustrated in Fig.2(b). Two pair of Weyl points are located exactly in the = 0 plane 
and another four pairs of Weyl points are located off the = 0 plane. Considering the 
four-fold rotational symmetry, it is then easy to understand that there are totally 12 pairs 
of Weyl points in the whole BZ. The appearance of Weyl points can also be derived from a 
k ■ p model with different types of mass terms induced by SOC, which will be introduced in 
detail in the Appendix. The band structures for the other three materials TaP, NbAs and 
NbP are very similar. The precise positions of the Weyl points for ah these materials are 
summarized in Table. I. 


FERMI ARCS AND SURFACE STATES 

Unique surface states with unconnected Fermi arcs can be found on the surface of a 
WSM. These can be understood in the following way: for any surface of a WSM, we can 
consider small cylinders in the momentum space parallel to the surface normal. In the 3D 
BZ, these cylinders will be cut by the zone boundary and their topology is equivalent to 
that of a closed torus rather than that of open cylinders. If a cylinder encloses a Weyl point, 
by Stokes theorem, the total integral of the Berry curvature (Chern number) of this closed 
torus must equal the total “monopole charge” carried by the Weyl point(s) enclosed inside. 
On the surface of the material, such a cylinder will be projected to a cycle surrounding 
the projection point of the Weyl point and a single Fermi surface cut stemming from the 
chiral edge model of the 2D manifold with Chern number 1 (or -1) must be found on that 













FIG. 3. Berry curvature from pairs of Weyl points, (a) The distribution of the Berry 
curvature for the kz — 0 plane, where the blue and red dots denote the Weyl points with chirality 
of +1 and -1, respectively; (b) same with (a) but for kz — 0.5927r plane. The insets show the 3D 
view of hedgehog like Berry curvature near the two selected Weyl points. 


circle. By varying the radius of the cylinder, it is easy to show that such FSs must start 
and end at the projection of two (or more) Weyl points with different “monopole charge” 

, i.e. they must be “Fermi arcs”. [6, 8, 15] In the TaAs materials family, on most of the 
common surfaces, multiple Weyl points will be projected on top of each other and we must 
generalize the above argument to multiple projection of Weyl points. It is easy to prove 
that the total number of surface modes at Fermi level crossing a closed circle in surface BZ 
must equal to the sum of the “monopole charge” of the Weyl points inside the 3D cylinder 
that projects to the given circle. Another fact controlling the behavior of the surface states 


9 



























is the MCN introduced in the previous discussion, which limits the number of FSs cutting 
certain projection lines of the mirror plane (when the corresponding mirror symmetries are 
still preserved on the surface). 

By using the Green’s function method [5] based on the tight-binding (TB) Hamiltonian 
generated by the previously obtained Wannier functions, we have computed the surface 
states for both (001) and (100) surfaces. They are plotted in Fig. 4 together with the FS 
plots. On the (001) surface, the crystal symmetry is reduced to C 2 v leading to different 
behavior for the surface bands around X and Y points respectively. Along the F-A or F-T 
lines, there are two FS cuts with the opposite Fermi velocity satisfying the constraint from 
the MCN for FZN plane. In addition to the MCN, the possible “connectivity pattern” of 
the Fermi arcs on the surface has to link different projection points of the Weyl nodes in 
a way that obeys the chirality condition discussed in the above paragraph. For the (001) 
surface of TaAs, the connectivity pattern of the Fermi arcs that satisfy all the conditions 
discussed is not unique. However, due to the fact that all the projective points on the (001) 
surface are generated either by a single Weyl point or by two Weyl points with the same 
chirality, the appearance of “Fermi arcs” on the (001) surface is guaranteed. The actual 
Fermi arc connectivity pattern for (001) surface is shown in Fig. 4(b), obtained by our ab- 
initio calculation on a non relaxed surface described by the TB model. Changes of surface 
potentials or the simple relaxation of the surface charge density might lead to transitions of 
those Fermi arc connectivity pattern and result in topological Fermi arc phase transitions 
on the surface. A very interesting point of the (001) surface states is the extremely long 
Fermi arcs which cross the zone boundary along the X to M line. Comparing with other 
proposed WSM materials, the Fermi arcs in TaAs families are much longer, which greatly 
facilitates their detection in experiments. 

Comparing to the (001) surface, the surface states in (100) surface of TaAs are much 
more complicated as shown in Fig. 4(c) and (d). The biggest difference between the (100) 
and (001) surfaces is that all the projected Weyl points on the (100) surface are formed 
by a pair of Weyl points with opposite chirality, which does not guarantee (but does not 
disallow) the existence of the Fermi arcs. The only constraint for the (100) surface states is 
the nonzero MCN of the FZN plane, which generates a pair of chiral modes along the FT 
line, the projection of the mirror plane, as illustrated in Fig. 4(d). 
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FIG. 4. Fermi arcs in the surface states, (a) The Surface states for (001) surface; (b) the 
corresponding fermi surfaces on (001) surfaces; (c) The Surface states for (100) surface. The dots 
illustrate the projective points of the bulk Weyl points on the surfaces, where the color represents 
the chirality of the Weyl points (blue for positive and green for negative), small dot represents 
the single projected Wely point, and large dot represents two Weyl points with same chirality 
projecting on top of each other, (d) the corresponding fermi surfaces on (100) surfaces. All the 
dots here are the projective points for a pair of Weyl points with opposite chirality. 


CONCLUSION 


In summary, a family of nonmagnetic WSM materials are proposed in the present paper. 
Each material in this family contains 12 pairs of Weyl points, which appear due to the lack 
of inversion center in the crystal structure and can be derived from the nonzero MCN for 
two of the mirror invariant planes in the BZ. The surface states of these materials form quite 
complicated patterns for the Fermi surfaces, which are determined by the both the chirality 
distribution of the Weyl points and the MCNs for the mirror invariant planes. 
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APPENDIX 


I. k-p Model of a Nodal Ring and the Appearance of Weyl Points Due to Spin-Orbit 
Coupling 


We perform &.k -p analysis in the vicinity of the E-point in the 3D BZ. We show that (i) 
a nodal ring (a closed line of band-touching points) is protected in the presence of a mirror 
reflection symmetry and SU(2) spin rotation symmetry (when spin-orbit coupling is absent) 
and (ii) when spin-orbital coupling (SOC) is turned on and all crystalline symmetries are 
preserved, the nodal ring may be partially gapped with several distinct possibilities: (i) Weyl 
nodes on the (OOl)-plane, (ii) Weyl nodes away from the (OOl)-plane and (hi) nodal rings 
on the (OlO)-plane and (iv) full gap. The realization of these different possibilities strongly 
depends on the specific form of SOC terms. 

The nodal ring around E-point can be modeled by a two-band k ■ p-theory, in the absence 
of SOC: 

Hoik) = d^ik)a„ ( 1 ) 

i=x,y,z 

where di(k) are real functions, and k = ikx,ky,kz) are three components of the momentum 
k relative to E-point along [100]-, [010]- and [001]-axes, respectively. In Eq.(l) we have 
ignored the kinetic term proportional to the identity matrix, as it is irrelevant in studying 
the band-touching. The mirror reflection symmetry, denoted by Mqio, is represented by 
M = Uz'- in the absence of SOC, we can choose a mirror symmetry that squares to unity. 
The form of the mirror operator is chosen such that the two bands have opposite mirror 
eigenvalues, an information obtained from the ab-initio calculation. The mirror reflection 
dictates that 


MHoikz,, ky, kz)M ^ = Hoik^;, -ky, kz), 

which translates into 


dx,yikxT ky, kz) dx^yikx, ky,kz), 

dzikx^ ky, kz) dzikx, ky, kz). 


( 2 ) 

( 3 ) 

( 4 ) 


Eq.(3) states that on the plane ky = 0, only dz is nonzero, and hence generically the equation 
dzikx, 0, kz) = 0 will have codimension of one, i.e., nodal line solution. Symmetry preserving 
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perturbations involve gradually changing the forms of dj’s without violating Eq.(3), so the 
nodal ring is robust against them. 

Another symmetry is present at the E-point: a twofold rotation C 2 about [001]-axis 
followed by time-reversal symmetry. This is because the rotation sends the S-point to its 
time-reversal partner, and a further time-reversal sends it back. This symmetry may be 

represented by C 2 T = UtK, where Ut is any symmetric and unitary matrix, and K is 

complex conjugation. Without SOC, the rotation about [OOlj-axis and the reflection about 
(OlO)-plane commute with each other, so we require 

[a,,KUT]=0, (5) 

and we may choose Ut = Oz- This symmetry places additional constraints on dj’s: 

ky, k^) ct^^Hq {kx^ ky, k2)o'zt ( 6 ) 


or 


dy^zikxi kyi dy^zi^kx^ ky, kz)^ 

dx ( kx , ky , k^ ) dx ( kx , ky , k^ ). 


( 7 ) 


Eqs.(3,7) determine the general form of our k ■ p-model. 

Now we consider adding spin-orbit coupling terms while respecting the symmetries at 
E-point. We first need to determine the matrix representations of the generators of the 
little group, i.e., Mqio and C 2 * T. Considering spin degrees of freedom, we know that (i) 
a mirror reflection consists of a spatial reflection and a twofold spin rotation about the 
axis perpendicular to the reflection plane, (ii) a twofold rotation involves a spatial twofold 
rotation and a twofold spin rotation about the same axis and (iii) time-reversal symmetry 
involves complex conjugation and a flipping of the spin. Eollowing these facts we obtain the 
matrix representations: 


M = iaz ® Sy, (8) 

C2T = K<yz ® Sx- 

Notice that now = — 1 and C^j, = 1, as needed. With spin degrees of freedom, each 
band in the previous spin-orbit coupling free model in Eq.[l] becomes two bands, and nodal 
ring becomes a four-band crossing. In the vicinity of the nodal ring, the addition of SOC is 


16 


equivalent to adding coupling between different spin components, i.e., ‘mass terms’, to the 
previous model. Here the name ‘mass term’ simply means that these terms are not required 
to vanish at the nodal ring by any symmetry. 

The symmetry of the nodal ring is just mirror reflection. The mass terms hence must 
commute with mirror symmetry, and a generic term on the ky = 0 plane is given by 


Hm = mi{k)sy + m2(k)azSy + m3{k)a^s^ + m4{k)a^s^ + m3{k)ayS^ + mQ{k)ayS^. (9) 

Note that these mass terms are in general k-dependent, as their values may change as k 
moves along the nodal ring, but the C 2 * T symmetry makes them satisfy (on the ky = 0 
plane) 


?7ri,2,4,6(/^x, 0, k^) = 0, -k^) = 0, k^), (10) 

^3,5(^0;) 0 , k ^;') ^ 3 ,5(^0:) 0 , fc^) 0 ) k ^^. 

A complete analysis of the band crossing in the presence of all six mass terms is unavailable as 
the analytic expressions for the dispersion are involved. However, one may see the qualitative 
role played by each mass term by analyzing them separately. From Eq.(lO), we see that 7713^5 
are odd under k^ —k^, while the others are even. This indicates that only mi_ 2 , 4 , 6 -fernis 
are responsible for band crossings appearing on the kz = 0 plane, while the band crossings 
away from that plane are attributed mainly to the presence of ms^s-terms. 

At ky = 0 plane, mi^ 2 -terms commute with Hq, so these terms, if of small strength, will 
split the doubly degenerate nodal ring into two singly degenerate rings, but not open gaps. 
The equations for the two new rings are given by 

dz{kz:,0,kz)±mi^2 = 0- ( 11 ) 

One should note that when mi-term (m 2 -term) is added, the two rings are the crossing 
between two bands with same (opposite) mirror eigenvalues. Therefore, the two rings from 
adding rrii term are purely accidental, and the two rings from adding m 2 term are protected 
by mirror symmetry. 

Next we discuss the effect of 7774^6 terms, which should in combination gives rise to the 
pair of Weyl nodes on the kz = 0-plane shown in the paper. The dispersion after adding 


17 


m 4,6 terms is 



where d? = + dy d^. With some straightforward algebraic work, it can be shown that 

the equation £'(k) = 0 (band-touching) is equivalent to the following three equations: 

4 = 4 = 0, (13) 

j2 2 2 

dy = rriQ — m^. 

When I me I > |m 4 |, these equations have at least one pair of solutions on 4 = 0-plane 
symmetric about ky = 0 with codimension zero: they are Weyl nodes on the 4 = 0 plane. 
In our simulation, we found only one pair of Weyl nodes appear on this plane, which can 
only be understood if m 4 e are k-dependent. The equations 4 = 0 and d^ = 0 determines a 
closed loop on the 4 = 0 plane. At the same time dy{kx, ky,0) = i-^/mg — ml has solutions 
that are symmetric about ky = 0 . Since 4 ( 4,05 0) = 0 , the solutions do not cross the 
4 = 0 line if are constants. Therefore, the solutions must be two lines, which make 
four crossings in total with the solution to 4 = 0. However, let us recall that all mass 
terms can also contain a linear function in 4 , so it is possible that uiq — vanishes for 
a particular 4- At that 4, the solution to 4 = 4^6 “ ^4 = 0 is satisfied at 4 = 0. If 
the point ( 4 , 0 , 0 ) is inside the loop that solves 4 ( 4 , 4 ’ 0 ) = 0 , then there must be an odd 
number of crossings of 4 = a/ mg — m\ and 4 = 0 . 

We can also understand the pairs of Weyl nodes that are away from kz = 0-plane. 
Consider a coexistence of both m 4 - and ms-terms. The dispersion is given by 

E(k) = (4 ± m4)2 -K (4 ± -F 4- (14) 

Solving A'(k) == 0 is equivalent to solving 


dz = 0, 

(15) 

4 = ±m4. 


dy = ims. 


The last two equations together give 


4 = ^/vm4^f\, 

(16) 

4 = \/ rn^X/v, 

(17) 
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where = ukykz, = Xk^ and dy = vky. 

When vni/^X < 0, there is no solution. When vru/^X < 0, there are four sets of solutions 
for {ky, kz). We can substitute them into = 0 to obtain the four Weyl nodes observed in 
our simulation. When 7714 is small, the four Weyl nodes are close to the crossing point of 
the nodal ring and the = 0 plane. 

We summarize the roles played by different mass terms: mi^ 2 -terms split the nodal ring 
into two non-degenerate rings. With mi- (m 2 -)term alone the ring is the crossing of two 
bands with same (opposite) mirror eigenvalues, ms^s-terms gap the nodal ring except at 
ky = kz = 0. m 4 -term alone or coexisting with ms-term fully gaps the ring. m 4 -term 
coexisting with ms produces four Weyl nodes away from the kz = 0 plane, me-term creates 
a pair of Weyl nodes on the = 0 plane, symmetric about ky = 0. 

II. Distribution of nodal lines and Weyl points 

A movie demonstrates the momentum space distribution of nodal lines calculated without 
spin-orbit coupling (SOC) and Weyl points with SOC. The movie can be download here: 
http://pan.baidu.com/s/IpJwTsmB. 
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